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Abstract 

A mathematical framework is proposed to predict the features of the (5 5 7) lath 
transformation in low-carbon steels based on energy minimisation. This theory gener¬ 
ates a one-parameter family of possible habit planes and a selection mechanism then 
identihes the (5 5 7) normals as those arising from a deformation with small atomic 
movement and maximal compatibility. While the calculations bear some resemblance 
to those of double shear theories, the assumptions and conclusions are different. Inter¬ 
estingly, the predicted microstructure morphology resembles that of plate martensite, 
in the sense that a type of twinning mechanism is involved. 
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1 Introduction 

The present article proposes a theory that predicts the formation of habit plane normals 
very close to (5 5 7), observed in steels with low carbon content (less than 0.4% |MW92j ). 
In fact, all the predicted habit plane normals are almost exactl3|^(2 2 3). Widely accepted 
models that result in (5 5 7) habit planes are double shear theories, e.g. |AB69[ lRC70j . 
and some of the most accurate explanations are due to the algorithm developed by Kelly 
|Kel92j . These can be seen as generalisations of the so-called phenomenological theory of 
martensite most notably developed by Wechsler, Liebermann & Read |WLI1,53] to explain 
the (3 1015) habit planes in plate martensite and Bowles & MacKenzie |BM54b] who 
applied their theory to explain the (2 5 9) and (2 2 5) habit planes also in plate martensite. 

*konstantinos.koumatosOgssi.infn.it 
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^E.g. one of the predicted normals is m = (0.4813,0.4956,0.7228) with z (m, (2 2 3)) = 0.7° and 
IIm- (2 2 3) lb = 0.012, where || • lb denotes the Euclidean distance. 
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A short-coming of single/double shear theories is the lack of a selection mechanism 
that picks the right lattice invariant shearing systems (see e.g. |Kel92l Table 1]), in turn 
leading to a large number of input parameters. To overcome this, one approach is to 
only allow shearing systems that arise from mechanical twinning, cf. Indeed, in 

the context of single shear theories |WLR53] and |BM54aj made this assumption and 
proposed that the martensite plates consist of a “stack of twin-related laths”. As pointed 
out in |KK12l p. 381f], TEM investigations in |KN61j showed that even though “under 
the optical microscope there is little sign that this is the ease, [...] when such steels are 
examined with the transmission electron microscope, arrays of very thin {11 2}m twins are 
indeed found” - marking a significant success from a theoretical prediction to an observed 
feature in plate martensite. 

Regarding (5 5 7) habit planes in lath martensite, it can be shown that for any reason¬ 
able choice of lattice parameters (see also Figure]^, a single shear theory with shearing 
systems arising from twinning in bcc crystals cannot give rise to them. However, in this 
paper we show that by introducing another level of twinning (“twins within twins”) we are 
not only able to explain {5 5 7} habit planes but also predict them by showing that it is 
the only possible family of habit planes that satisfies a condition of maximal compatibility 
and a condition of small overall atomic movement. Under this interpretation each lath 
may be seen as a region of twins within twins. In other materials, twins within twins have 
commonly been observed purely in martensite |Arl90l IDCV’*'!^ as well as along interfaces 
with austenite [BKS10| . Moreover for lath martensite it has been observed in [SW83) 
that “Twinning within a lath may he heavy [...]. In any event, whenever an exact twin 
relationship was identified, it was found to he a result of twinning within a given lath and 
not of a twin relation between adjaeent laths. [...] It is believed that the existenee of heavily 
twinned loeal regions of laths, which may appear as separate laths in contrast images, may 
have caused some misinterpretation in earlier work on lath martensite. ” 

As in single shear theories, twinning of twins is macroscopically equivalent to a simple 
shear of a simply twinned system (cf. text below ([^) and thus the step from twins within 
twins is in analogy with the step from a single to a double shear theory. 

The strength of the theory presented here is that it enables one to predict (5 5 7) 
habit planes only assuming the lattice parameters of austenite and martensite. This is 
particularly striking when compared to the double shear theory in |Ke]92] where the 
“caleulation strategy was to select one of the possible S 2 [second shear] systems and then 
perform calculations for the Si [first shear] systems [...] over a range of values 0/52 [the 
shearing magnitude of S 2 ]. [...] The sign of g 2 was selected by trial and error depending 
on whether the habit plane moved towards or away from (5 5 7).” 

2 A model of phase-transformations based on nonlinear elas¬ 
ticity 

The theory proposed in this article is derived from the Ball-James model |B.I87j - based 
on nonlinear elasticity and energy minimisation - and expands on previous work by Ball &: 
Carstensen |BC97) on the possibility of nonclassical austenite-martensite interfaces. Even 
though lath martensite is commonly associated with a high dislocation density, slip and 
plasticity the present model does not take such effects into account (see also Section]^. 

Interestingly, the Ball-James model recovers the results of the phenomenological theory 
of martensite, as can be seen through a comparison of the derived formulae for the habit 
planes between twinned martensite and austenite (cf. |BJ87[ eq. (5.89)] and |WLR53l eq. 
(33)-(34)]). For a self-contained account of the phenomenological theory or of the Ball- 
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James model the reader is referred to the monoeraohs |Bha m and |Bha03] . respectively, 
both addressed to non-specialists. 

In the Ball-James model, which neglects interfacial energy, microstructures are iden¬ 
tified through minimising sequences y^, k = 1 , 2 ,..., for a total free energy of the form 

E{y)-.= j^We{Dy{x))dx. (E) 

Here, is a region representing the reference configuration of undistorted austenite at 
the transformation temperature and y{x) denotes the deformed position of particle x in 
n. We remark that passing to the limit in these minimising sequences, corresponds in a 
very precise way to passing from a micro- to a macroscale, so that the limits themselves 
can be identified with the macroscopic deformations. The energy density W 0 {F) depends 
only on the deformation gradient F = Dy, a 3x3 matrix with positive determinant, and 
the temperature 6 . Also Wq is assumed frame indifferent, i.e. Wq^RF) = W 0 {F) for all 
rotations R - that is, for all 3x3 matrices in SO(3) = {R ■ R^R - I, deti? = 1} and must 
respect the symmetry of the austenite, i.e. W 0 {FP) = W 0 {F) for all rotations P leaving 
the austenite lattice invariant. For cubic austenite there are precisely 24 such rotations. 
Below the transformation temperature, Wq is minimised on the set K of martensitic energy 
wells, that is W 0 {F) is minimal for E e A' = uili SO(3)t/i. The 3x3, positive-definite, 
symmetric matrices Ui are the pure stretch components of the transformation strains 
mapping the parent to the product lattice. For example, in the case of fee to bcc or fee to 
bet, these are given by the three Bain strains 

Ui = Bi = diag(? 72 ,r?i,r/i), U 2 = B 2 = diag(r?i, 772 , m), U 3 = B 3 = diag(r7i, 771 , 772 ), 

where 771 = and 772 = ^. Here oq is the lattice parameter of the fee austenite and a, 
c are the lattice parameters of the bet martensite {a = c for bcc). The notation Bi, B 2 , 
i ?3 has been chosen to emphasise that we are in the Bain setting and to stay consistent 
with the literature. We remark that the Bain transformation [BD24] is widely accepted 
as the transformation from fee to bet/bcc requiring least atomic movement; for a rigorous 
justification see [KM] . 

A convenient way to understand the relation between microstructures and minimising 
sequences is illustrated by the following example (cf. |BJ87] 1. 

Example 1. (Austenite-twinned martensite interface) 

Suppose that a region of martensite is occupied by an array of twin related variants 



Figure 1: An energy minimising sequence modelling twinned martensite. The limit k 00 
corresponds to an IPS leaving the plane with normal m invariant. 

Ai and A 2 with relative volume fractions 1 - A and A. The strains Ai and A 2 in general 
cannot be invariant plane strains (IPS), equivalently, they cannot form a fully coherent 
interface with austenite, represented in this model by the identity matrix I. However, for 
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specific volume fractions (given by for the Bain strain), the average deformation 
strain of the twinned region (1 - A*)^i + X*A 2 may indeed become an IPS. In terms of the 
nonlinear elasticity model, this inability to form a fully coherent interface at the micro¬ 
scopic level, implies that the austenite-twinned martensite configuration cannot exactly 
minimise the energy (|^. Nevertheless, one can construct an energy minimising sequence 
y^, k = 1,2,..., with gradients Dy^ as in Fignre[^ The limit of this sequence is precisely 
the average strain/total shape deformation (1 - A*)^i + \*A 2 that is an IPS. Although 
this average strain does not minimise the energy (j^, it can be interpreted as a minimiser 
of a corresponding macroscopic energy. At the microscopic level, one would observe some 
specific element of the above minimising sequence rather than the limit, due to having 
neglected interfacial energy. 

The set of all matrices that can arise as limits of minimising seqnences for the energy 
( [^ is referred to as the quasiconvex hull of the martensitic energy wells K, denoted by 
Hence, the set corresponds to all possible homogeneous total shape deformations 
that are energy minimising at the macroscopic level. Then, the reqnirement that (np to an 
overall rotation R) a martensitic microstructure with a total shape deformation F e 
is a strain leaving the plane with normal m invariant, amounts to finding a vector b such 
that 

iiF = I + b0me (1) 

Here, b 0 m denotes the 3x3 matrix (b 0 m)jj = binij. Writing b = b/|b|, yields the 
eqnivalent expression RF = I + |b|b 0 m, implying that RF is the IPS given by a shear 
on the plane with normal m, with shearing direction b and shearing magnitude |b|. In 
particnlar, if the transformation from parent to product phase is volume-preserving, F is 
a simple shear, corresponding to the vectors b and m being perpendicnlar. We note that 
Q equivalently says that RF can form a fully coherent planar interface with austenite of 
normal m. Also by frame-indifference, the austenite can be represented by any rotation. 
Hence, a further rotation Q of the martensite resnlts in QRF = Q + (Qb) 0 m, so that 
the rotated martensite can still form a fully coherent interface with austenite of the same 
normal m. That is, rotations from the left cannot change the habit plane normal. 

3 Comparison to the phenomenological theory 

A common feature in both the phenomenological theory and the Ball-James model is to 
construct (up to an overall rotation R) a total shape deformation F that is an IPS. In the 
literatnre, various algorithms have been proposed for the calculation of the corresponding 
elements of the shear, i.e. the magnitnde, direction and normal (cf. Q). For example, 
see |MB54l IWLR53] in the context of twinning/single shear theories and |AB69l IRC 70 j 
for double shear theories. 

In general, the problem of finding an overall rotation R and shearing elements snch 
that RF = I -H b 0 m can be simplified by only considering the Cauchy-Green strain tensor 
C = {RF)^{RF) = F"^F and thus factoring out the overall rotation R. The following 
Proposition ( |B.I87t Proposition 4]), allows one to calculate the shearing elements b and 
m in terms of the principal stretches and stretch vectors of F. The overall rotation R can 
then be found by substitnting b and m back into equation Q. 

Proposition 1. Let C be a symmetric 3x3 matrix with ordered eigenvalues Ai < A 2 < 
A 3 . Then C can he written as 


C = (I -H m 0 b)(I -H b 0 m) 
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for some b, m if and only if Xi > 0 and A 2 = 1. Then, there are at most two solutions 
given by 

b = x-1 , 

m = p~^ kv^As - IV 3 ), 

where p 0 is a normalisation eonstant, k e {- 1 , 1 } and vi,V 3 are the (normalised) 
eigenveetors of C corresponding to Ai and A3. 

A short interlude on martensite twins 

In the material science literature twins are often described as two phases related by a 
specific 180° degree rotation or, equivalently, a reflection. In the mathematical literature 
a twin is usually characterised by the existence of a rank-one connection between the two 
deformation strains Ai, A 2 corresponding to the two phases, i.e. the existence of vectors 
a and n such that 

A 2 = Ai + a 0 n. 

A fully coherent interface between the two phases is then given by the plane of normal 
n. This is because for any vector v on that plane, i.e. v • n = 0, we obtain A 2 V = 
AiV + (v • n)a = Ai v. Also, note that A 2 = (I + a 0 A}^n) Ai so that the lattice on the one 
side of the interface can be obtained by shearing the lattice on the other side along the 
twin plane , in the shearing direction with shearing magnitude |a||A{^n|. The 

latter expression enables one to calculate the vectors a and n by Proposition through 
the identification F = A 2 A{^, that is the relative deformation between the two phases is 
an IPS. In view of single shear theories, the above expression can equivalently be written 
as A 2 = Ai(I + A{^a 0 n) and thus A 2 can be obtained as a shear of the parent lattice, 
followed by Ai. 

Hence, in the case of twins between two martensitic energy wells SO{3)Ui and SO{3)Uj 
one needs to solve the equation 


QUj = [/j + a 0 n 


( 2 ) 


for the rotation matrix Q and the twinning elements a and n. If the transformation 
strains Ui and Uj are related by a 180° rotation, this calculation simplifies significantly by 
Mallard’s Law (see [Bha03l Result 5.2] or below). In particular, this assumption holds for 
Ui = Bi and Uj = Bj, i.e. for the Bain transformation from fee to bct/bcc. 

Proposition 2. (Mallard’s Law) 

Let U and V satisfy V = PUP for some 180° rotation P about a unit vector e, i.e. 
P = -I + 2e 0 e. Then the equation QV = t/ + a 0 n admits two solutions given by 


a= 2 


U-^t 


\u 


-To|2 


Ue\, 


2N 


;Ue, 


|[/e|2 

In each case, Q = (1/ + a 0 n)P“^. 


n = e. 


(I) 


n= ^{\Ue\^e-U^Ue). 


( 11 ) 
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We conclude this interlude by remarking that twins described by the first solution in 
Mallard’s law are Type I twins and the corresponding lattices are related by a 180° rotation 
about the twin plane The second solution in Mallard’s law describes Type II twins 

and the lattices are related by a 180° rotation about the shearing direction |^. It may 
happen, and it does for the Bain strains, that there are two rotations by 180° relating 
U and V. In this case, there are seemingly four solutions from Mallard’s Law, however. 
Propositionsays that there cannot be more than two. Indeed, the Type I solution using 
one 180° rotation is the same as the Type II solution using the other 180° rotation and 
vice versa. In particular, the lattices on either side of the interface are related by both 
a 180° rotation about the twin plane and a 180° rotation about the shearing direction. 
Solutions of this type are compound twins. 

3.1 Single shear theories 

Henceforth, we only consider the Bain strains Bi, B 2 and B 3 for the fee to bct/bcc 
transformation in steel. In single shear theories the total shape deformation is assumed 
to be decomposable into F = RBS where i? is a rotation, B is one of the Bain strains and 
S = I+d0p is a shear whose specific form varies in the literature. In the Ball-James theory, 
the total shape deformation F must be macroscopically energy minimising, thus restricting 
the form of the shear S. The most important case is when S arises from twinning. As 
in Example the average strain corresponding to a twinning system between Ai = Bi 
and A 2 = QB 2 , satisfying with volume fractions 1 - A and A, respectively, is given by 
(1 - \)Bi + XQB 2 = Hi + Aa 0 n, for any A e (0,1), where the elements Q, a and n can be 
calculated by Mallard’s Law (cf. Proposition]^ applied to U = Bi with either e = (1,1,0) 
or e = (1, -1,0), i.e. the resulting twins are compound. By simple algebraic manipulation, 
the average strain can be written as 

Bi(I + ABi^a0n) = BiAa, (3) 

i.e. a single shear Sx = 1 + XB^^a 0 n, with det Aa = 1 , of the parent lattice followed by the 
Bain strain Bi. 

By Proposition]^ to make the total shape deformation Fx = RBiSx an IPS, the volume 
fraction A needs to be chosen such that the middle eigenvalue of F^Fx is equal to one. 
In particular, since one of the eigenvalues must be made equal to one, the expression 
det{F^Fx - I) must vanish, giving rise to the two solutions A* and 1 - A*, where 

= 2 ^ 2 '\/(‘^-'d2-vl){vl-‘^vh2+vl)- (4) 

It is important to then check that it is indeed the middle eigenvalue that is equal to one. 
For each of the values A*, 1 - A*, we can calculate two habit plane normals according to 
Proposition]^ One of these normals is, up to normalisation, given by (hkl), where 

By considering the remaining normals and all possible pairs of twin related Bain strains, 
one recovers the entire family of normals {hkl}. In Figure]^ the components of one 
of these normals are plotted for a typical range of lattice parameters r/i, r /2 in the fee to 
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bct/bcc transformation. We immediately note that for rji = 1.1 and 772 = 0.86 the predicted 
habit plane normal arising from simple twinning is almost exactly (3 10 15) and for rji = 
1.11 and r ]2 - 0.86 the habit plane normal is almost exactly (2 5 9). The corresponding 
ratios of tetragonality are given by c/a = ^ 2771/772 ~ 1-105 and rj 1.095 respectively. This 
in excellent agreement with the observations in e.g. |RR56j of (2 5 9) habit planes in steel 
with carbon content in the range 1.4-1.8 wt-% as well as the theoretical and experimental 
results in e.g. |KN61j . |WW71| and |GT49j of (31015) habit planes in highly tetragonal 
martensite. 



Figure 2: Coordinates of habit plane normal (blue, red, green) arising from simple twinning 
for different values of 771 , 772 . 


3.2 Double shear theories 

Similarly, in double shear theories the total shape deformation is assumed to be decom¬ 
posable into F - RBS 2 S 1 where R is a rotation, B is one of the Bain strains and ^i, 
S 2 are two shears. Above we have seen how twinning can be regarded as an instance of 
a single shear theory. In this section we show how an additional level of twinning, i.e. 



Figure 3: A configuration of twins within twins macroscopically leaving the plane with 
normal k invariant. 

twins within twins, results in an instance of a double shear theory, consistent with the 
Ball-James model. To visualise this type of microstructure, we revisit the construction in 
Example to construct two twinning systems one with Ai = Bi, A 2 = QB 2 and average 
strain Ri - 1 - Aa 0 n and another one with say Ai = Bi, A 2 = Q'B^ and average strain 
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Bi + Aa' 0 n'. For each A there exists a rotation R'^ such that 
R'y^{Bi + Xa' 0 n^) = {Bi + Aa 0 n) + 0 


and thus the two twinning systems are macroscopically compatible with a fully coherent 
interface of normal between them (cf. Figure^. The elements R'^, b;^ and can 
be calculated by Mallard’s Law (cf. Proposition!^ applied to U = i?i + Aa 0 n with 
e = (0,1,1), giving rise to a Type I and a Type II solution. Unlike in the single shear 
theory, these twins are not compound and therefore we distinguish these two solutions by 
the superscript a, a = 1,2 for Type I and Type II respectively. Finally, it can be shown 
that the volume fractions A in each of the twinned regions must necessarily coincide. As 
in Example we can construct an array of twins between the two twinned regions with 
respective volume fractions 1 - /r and /x and average strain given by 


Si + Aa 0 n + |ub“ 0 m", A,/xe(0,l). 


(5) 


Simple algebraic manipulation allows one to write ([^ as SiS^(A,/u)Si(A) where Si = 
S'i(A) = I + ASj(^a 0 n and S 2 = (A, ^) = I + /xSj'^b" 0 and thus an instance of a 

double shear theory. We note that detSi = detS 2 = 1. 

By Proposition!^ in order to make the total shape deformation = SSiS"(A, /x)5i(A) 
an IPS, the volume fractions A, /x need to be chosen such that the middle eigenvalue of 


is equal to one. Solving 


det(F,“ 


OL T ptCt 


•I) = 0 


for each fixed A, gives rise to a quadratic equation in /x for each choice of a, which can 
be solved explicitly (see for the full details). The expressions are lengthy and we 

refer to Figure which visualises the dependence /x“(A) for the volume-preserving Bain 
strain. The endpoints of these curves correspond to the vanishing of one of the twinning 
systems (cf. (|^) and hence to the collapse of the system of twins within twins, to a simple 
twinning system with volume fractions given by Q- The figure remains qualitatively 
the same for typical lattice parameters. For each a and each admissible pair (A,;u“(A)) 
with corresponding strain we can calculate two one-parameter families of habit 

plane normals through Proposition!^ By considering all possible combinations of twinning 
systems in our construction, we then obtain all crystallographically equivalent normals. 
For the volume-preserving Bain strain, the habit plane normals that can arise from the 


family F) 


A,/ii(A) 


are visualised in Figure 


However, due to algebraic complexity, it is 
difficult to write down a formula for the habit plane normals with an explicit dependence 
on r/i, r \2 and A. 


Why twins within twins? 

It is natural to assume that the observed total shape deformation F requires small overall 
atomic movement (see also |KMj ) relative to the parent phase of austenite. A measure of 
this distance is the strain energ-^ given by 

dM(F,I) = |F^F - Ip = - I)^ 

i=l 

where |Ap = Tr(A^A) denotes the Frobenius norm and i'i{F) the principal stretches of F. 
It can be shown that any microstructure with small strain energy must necessarily involve 


^Alternatively, one may use dM(F, I) = “ 1)^ which yields the same results. 
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Figure 4: Ternary plot of volume fractions 1 - A of i?i, (1 - ^“(A))A of B 2 and /r"(A)A of 
i ?3 that make the twins within twins an IPS for 771 = 2^^® and 772 = 


all three Bain variants in roughly similar volume fractions. In particular, this cannot be 
the case for an array of twin related variants and we ought to consider at least twins 
within twins (see also Fig. [^. Although, introducing even further levels of twinning can 
reduce the strain energy, one could argue that interfacial energy contributions, which are 
not accounted for in this model, may inhibit such behaviour. 


4 A new theory for the (5 5 7) lath transformation 


Combining the fact that {5 5 7} cannot result from simple twinning (cf. Fig. and 
that twins within twins are preferable in terms of strain energy, we build a theory that 
predicts {5 5 7} habit plane normals solely based on energy minimisation and geometric 
compatibility. 

Firstly, the one-parameter families of habit plane normals obtained from twins within 
twins (see Section 3.2), contain normals very close to any {5 5 7}. This is at least the 
case for lattice parameters close to 771 = 2^/® r; 1.12, 772 = r; 0.79 corresponding to a 
volume-preserving transformation from fee to bcc. This regime of parameters is suitable 
since {5 5 7} habit planes are observed in low-carbon steels where the transformation is 
very nearly fee to bcc. The resulting one-parameter families of habit plane normals, along 
with their crystallographically equivalent ones, are shown in Figure We stress that the 
only free parameter in the generation of these normals is A which fixes the choice of the 
shearing systems, based only on the energy minimising property of the microstructure. 
Secondly, out of these one-parameter families of normals, our theory can identify the 
{5 5 7} habit plane normals as those satisfying a criterion of maximal compatibility. To 
this end, revisiting our construction of twins within twins there is a choice (cf. ([^) of using 
either average strain, corresponding to the Type I solution from Mallard’s 

Law, or F^^2^x) corresponding to Type II. Figure ^ shows the strain energy associated 
with the two macroscopic strains as a function of A. It is clear that the strain energy 
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Figure 5: Possible habit plane normals for and r /2 = 

2-1/3 _ Yellow points correspond to habit plane normals arising from simple twinning. 


of -Fx^i(_x) i® significantly smaller than that of F‘^^2 (^x) is thus preferable. We also 
note that, in agreement with the previous section, the strain energies of both and F'^ 
increase rapidly as the volume fraction of Bi approaches 0, A* or 1 - A*, that is as the 
microstructure reduces to a single twinning system. Further, we remark that the ^2(^x) 
with minimal strain energy result in habit plane normals which are very nearly {111} 
(see also Fig. 5). Nevertheless, the strain energies of any of the F(x/ii(A) that give rise 

Even though, the strains resulting in {5 5 7} 


to {5 5 7} normals are lower (cf. Fig. [^. 



volume fraction oi Bi 


Figure 6: Strain energies for r/i = 2^^® and r }2 = 2 


habit plane normals do not minimise the strain energy, they satisfy a strong criterion of 
compatibility. To understand this one must think in terms of the dynamic process of 
nucleation. As austenite is rapidly quenched, the martensite phase nucleates at various 
sites. The strain in a given nucleation site may need to be an IPS but otherwise, has no 
reason to be the same as the strain in any other site. Nevertheless, there are essentially 
only three distinct families of systems of twins within twins and these can be classified by 
the Bain variant which is present in both of the simple twinning systems that comprise 
the overall microstructure. In Figure these three families are distinguished by colour. 
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As the nuclei grow and approach other nuclei, they need to remain compatible with each 
other. Remarkably, the only habit plane normals that arise from deformations with low 
strain energy and can be reached by all three families are {5 5 7}. In Figure]^ this can be 
seen from the fact that all differently coloured curves intersect close to (5 5 7). For any 
two such regions of twins within twins with corresponding average strains I + bi 0 (557) 
and I + b 2 0 (5 5 7), one can see that 

(I + bi 0(557))-(I + b2 0 (5 57)) = (bi -b2)0(557), 

implying that they can meet along a fully coherent planar interface of normal (5 5 7). Of 
course, any nucleus interacts with its neighbours faster than it does with distant nuclei. 
As a result, blocks of similarly oriented regions of twins within twins (laths) may form 
whose overall orientation may differ from that of other blocks. 



Figure 7: Coordinates of the habit plane normals (blue, red, green) for the macroscopic 
strains F'x ^ ^ (^f- Proposition 1). The lattice parameters are rji = 2^/® 

and r ]2 = 2'^/^, corresponding to the volume-preserving Bain strain. 


5 Concluding remarks 


The theory proposed here for the prediction of (5 5 7) habit plane normals has two possi¬ 
ble interpretations. On the one hand, it can be seen as a purely macroscopic theory. In 
particular, it is an instance of a double shear theory with a precise algorithm to produce 
the required shears based on energy minimisation, without the need for any further as¬ 
sumptions. All possible habit plane normals that can arise from the one parameter family 
(indexed by A) of macroscopic deformations are shown in Figure Table jl then 

lists the elements of the twinning systems for the values of A that produce a near (5 5 7) 
habit plane. With the help of ([^ it is easy to convert between the twinning and shearing 
systems and thus compute the elements Si and S 2 required in a double shear theory. At 
this macroscopic level it is not possible to distinguish between twins within twins, a single 
twin and one slip system, and a single variant and two slip systems. 

On the other hand, a physical mechanism for the formation of (5 5 7) habit plane nor¬ 
mals is proposed and thus a specific morphology on a microscopic level. According to this 
interpretation, each lath may itself be a region of twins within twins with a corresponding 
lath boundary of normal (5 5 7). This type of morphology is depicted in Figure]^ with k 
being a {5 5 7} normal and the other elements are as in Table This morphology is a 
direct consequence of the underlying theory and it would be very interesting if it could be 
put to experimental scrutiny. 
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A 

0.576 

0.659 

0.762 

a 

[.374,-.529,0] 

[.374,-.529,0] 

[.374,-.529,0] 

n 

2-i(l,l,0) 

24(1,1,0) 

24(1,1,0) 

p\X) 

0.581 

0.621 

0.546 

bA 

[.130, .234, .315] 

[.135, .260, .359] 

[.137, .289, .412] 

m 

24(0,1,-!) 

24(0,1,-!) 

24(0,1,-!) 


Table 1: Elements of the twinning system leading to {55 7} habit plane normals. 
The remaining {5 5 7} normals can be obtained from the crystallographically equivalent 
systems. 
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